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1 INTRODUCTION 



This paper addresses the question of lattice sizes needed in practical applica- 
tions in order to fulfil the asymptotic scaling behaviour predicted by Borgs 
and Kotecky 0, close to a first order phase transition point. 

• Very strong first order transitions are easy to detect. The system be- 
havior is very close to the infinite volume behavior. Ergodicity is broken. 
Thermo dynamical quantities are discontinuous at the transition point, with 
metastable branches. A starting configuration half ordered, half disordered 
will relax to very different states on both sides of the transition. 

• In less clear cases, one must simulate systems of increasing volumes 
L'^ and try to convince oneself that the above described very large volume 
behavior is approached. Let us use the language of energy driven transitions 
in what follows, and introduce the energy probability distribution Pl{E). In 
the transition region, it has two peaks of heights P£ (ordered phase) and 
Pi (disordered phase), separated by a minimum of height P|"*". For a first 
order transition, at fixed P£/P^, p™"/p£ decreases as L grows. This has 
been proposed by Lee and Kosterlitz as an indicator of first order phase 
transitions (for a microcanical version see 0, 0). One can also plainly look 
at a plot of Pl{E) and decide "by eye" whether it looks more and more like 
two delta functions as L grows. Another class of indicators are moments 
of the distribution Pl{E), that go to zero in the large volume limit, for all 
temperatures, but at a first order transition point. An example is the energy 
fluctuation CV/L'^ = l3'^{< E"^ > — < E >^), another is Binder's cumulant 
BL = \{1- <E^> / <E^ >2). 

• The most sophisticated (and trustworthy) method is finite size scaling. 
One insists in seeing the finite size behavior as predicted by the theory in the 
vicinity of a first order phase transition. In that case only one can be pretty 
sure that the trend observed for lattices of increasing size does continue up 
to the thermodynamical limit. One insists on seeing 

1 TDmin ( T \ 

l.lnf--^r^A/L (1) 

BLmm--BL^^^ ^BL^'^^IL'' (3) 
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where CVmax is the maximum of the specific heat, and BLmin is the minimum 
of BL. 

The above described large volume behavior has never really been observed 
in a numerical simulation. The results of all high precision simulations show 
strong curvatures in e.g the curve representing BLmin as a function of 1/L'^. 
This is the case of the 3-d 2'(3) Potts model with antiferromagnetic admixture 
on lattices up to 48^ 0. This is also the case of the 5 states 2-d Potts model 
[0] on lattices up to 256^. In the later case, the exact asymptotic limit 
is BLmin ~ —.44 10^'^ whereas the data presented would suggest a much 
lower value (~ —1.2 10^'^). This shows that the BLmin indicator is not a 
trustworthy indicator of first order phase transitions on lattices up to 256^ 
when the correlation length is equal to a few thousands lattice spacings. 

Results of a simulation of the g = 10 2-d Potts model on lattices as large 
as 50^ can be found in Ref. This is a model with a strong, obvious, first 
order transition, with an exactly known correlation lengthQ^(/3f") = 10.56. 
The results for the extrema of CV/L'^, BL and f/4 || 

UA = <iE-<E>r> (4) 
<{E-<E >)2 >2 

and for the value CV{l3t) / L'^ are compared with the large volume predictions, 
namely the extrema behave like X^^) + X^^) /L"^ + ^(l/L^'^), and CV{f3t)/L'^ 
behaves like X^^) + X^^)/^'^ + 0{e'^^) for some b > 0. The four different 
constants {X*^^)} are exactly known. One single unknown parameter, e.g. the 
ordered specific heat at the transition point Co, determines the four X^'^^^s. 
The precision of the data is very good, and in all cases strong deviations 
from the X^^-* +X^'^^ /L'^ limiting behavior are seen. These corrections do not 
seem to behave simply as function of L, and are definitely not under control. 
Not surprisingly, the values for the four slopes {X'^2)j one would infer from 
the data give inconsistent estimates^ of Co- Furthermore, the data do not 
substantiate the prediction that CV{Pt) is asymptotic earlier than CVmax- 

^The correlation length of the 2-d Potts model at the transition point can be computed 
exactly using the 6- vertex representation. The result is interpreted as the disordered phase 
correlation length in the limit /? ^ (3^ . The ordered phase correlation length is smaller. 
Duality does not relate theses two correlation lengths. 

^Due to a systematic typing error, every times a numerical value for Cd is given in Ref. 
Q (e.g. in Table 1), the value for the ordered phase Co is meant. In particular, the q = 10 
curves are drawn using the estimate of Co = 12.7 ± .3, not Cd = 12.7 ± .3. 
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In conclusion, the FSS behavior is not observed even for a transition as 
strong as the 2-d g = 10 Potts model transition, on lattices as large as 50^. 
Three possible explanations come to the mind: i) Eqs.0,|] are only proven 
in the large q limit, they may not hold down to g = 10. ii) Much larger 
lattices may be needed in order to extract the true asymptotic behavior, 
even though on the largest lattices simulated Pl{E) has a textbook first or- 
der shape, iii) A programming error may be possible. We are however in 
the position to compare high precision results for the q = 10 Potts model, 
which were obtained in two completely independent simulations |Q |10|, e.g. 
using different random number sequences. Results for thermodynamic quan- 
tities agree on same lattice sizes to high precision, making the existence of a 
programming error unlikely. 



2 Exact Results 

Those have been obtained |jl| for models that can be represented by a contour 
expansion with small activities, like the q states Potts model for large q. 
In such a case, the partition function for a L'^ lattice with periodic boundary 
conditions can be written as 

Z{P, L) = e-^'P^-^^P) + qe-^'i^UP) + 0(e-''^)e-^^(^)^' ; 6 > (5) 

where /o(/5) and are smooth L independent functions. The free en- 

ergy is /(/?) = min{/o(/3), /<i(/3)}. The phenomenological two gaussian peak 
model of the energy probability distribution Pl{E) introduced by Binder and 



Landau |]TT], [T^ follows through inverse Laplace transform. The above exact 
result fixes the relative weights of the two peaks: At the infinite volume limit 
transition point, P = Pt, the ordered and disordered peak weights are exactly 
in the ratio q to one. 

To the order in l/L*^ which we consider, all quantities are expressed in 
terms of Pt and of the energies and specific heats of the two coexisting phases. 
The transition temperature. Eg, and the difference Co — Cd are known 



exactly for the 2-d Potts models |13]. Values for the g = 10 and g = 20 cases 



can be found in Tab.|^. It follows from Eq.^ that the specific heat 

CV = I3''L\< E^>-<E>^) (6) 
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has a maximum at 

5(2) 
'cv 

The height of this maximum increases hnearly with L'^ 



^{CV^a.) = A - + % + 0{l/L^% (7) 



CV^a. = L'^iE, - Eaf + CV^(2) + 0{l/L'). (8) 

whereas for fixed /5 7^ CV{[3) goes to a constant, as L goes to infinity. 
One finds that BL reaches a minimum equal to |I|, |15| 

at the point 

= fl - '"''ifl^f f + ^ + 0( V^"). (10) 

Expressions of the coefficients BL^'^\ and CV^^-* as functions of 
the Ei^s and Q's can be found in [|l^]. f/4 reaches a minimum p| 



at the point 



/3(f/4_)-A-;^— ^;^+ ^2.2/3|(E„ - E,)3 +^^^/^ )• ^^2) 

The above formulae for the extrema CVmax, BLmin, UAmin and the cor- 
responding effective /3's have higher power law corrections that may hide the 
asymptotic behavior on lattices that can be simulated. In contrast, the ex- 
pressions for bulk averages evaluated at the (infinite volume limit) transition 
point /3 = /St do not have power law corrections, as a consequence of Eq.|], 
e.g. the specific heat behaves like 

CV{P,) = + - Edf^^ + ^(^"'") (13) 
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3 Simulation of the q=20 model using the 
Multicanonical Algorithm 

We have decided to simulate the q = 20 Potts model, as an example of 
a model with much stronger phase transition than the g = 10 model (see 
Tab.l). The conventional Metropolis (and Swendsen-Wang |TB]) algorithm 
suffers from exponential slowing down close to a strong first order point, 
namely the probability to jump from one phase to another goes like Po.d ^ 
g-2X,di'' ^ ^ where A"'^ is the order-disorder surface tension. This makes sim- 
ulation of a ~ 50^ lattice with such algorithms impossible even with vastly 
more powerful computers than available now. 



q 


Co 


Ed — Eo 


Cd — Co 




{Co/{Eo-Edff'' 


10 


12.7 ± .3 


0.69605 


.44763 


10.56 


5.12 ±.06 


20 


5.2 ± .2 


1.19416 


.77139 


2.70 


1.91 ± .04 



Table 1: The g = 10 and q = 20 2-d Potts model: Estimated values of 
Co, exact values for Ed — Eo, Cd — Co, and correlation length ^(/3 Pf~). 
A transition is strongly first order, on L'^ lattice when L » C,, L'^ » 

Col {Ed — EoY 

Such a simulation became possible with the invention of the "multicanon- 



ical ensemble" |jTO[ . It amounts to perform the simulation with a weight fac- 
tors designed in such a way that Pl{E) is smooth between Eo and Ed, and to 
"reweight" the events when computing expectation values. The multicanon- 
ical algorithm has only polynomial slowing down. 

The multicanonical algorithm is only partially vectorizable. On a 38^ 
lattice, our program takes 7.7 micro-sec to update a spin on a CRAY X-MP, 
27.8 micro-sec on a RS6000-550 and 35.0 micro-sec on a DEC 6000-510. We 
have made runs on lattices indicated in Tab.|^. Most runs were performed on 
workstations, using the pseudo random number sequence Jj = li-uj ® Ii-250, 
where ® is the logical exclusive OR. Run ^4 at L = 24, and run #1 at 
L = 38 were performed on a CRAY X-MP, using the congruential pseudo 
random number generator RANF(). Statistical analysis was performed using 
Jackknife technique, with first order bias correction. In many cases, we 
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performed several simulations for a given lattice size, using more or less 
optimized weights. Each run gave us results for e.g. CVmax- Those were 
averaged using weighted averages. In some cases, we obtained ^ bad as 
~ 2 p.d.f. We believe this would disappear with higher statistics. 

Our results for CVmax/ L'^, BLmm and Uimin can be found in Figs.|l], || 
and together with the theoretical estimate using the value Co = 5.2 ± .2. 
In all cases the predicted X^^^ + X^^^/L'^ behavior is approached on the 
largest lattices, furthermore the values of the three slopes {X*^^)} one infer 
from our data give consistent estimates of Co- This was not true in case 
of the g = 10 simulation on similar size lattices. Our results for P (CVmax) 
and f3{BLmin) can be found in Fig.^. Perfect agreement is found with the 
theoretical predictions. The effect of the correction is not visible, and 

the data for ^{Uimin) (not plotted) are not distinguishable from those for 
P{CVmax)- The behavior of CV{Pt) is shown in Fig.^ No deviation is seen 
from the predicted X^^^ + X^^y behavior. This agrees with the prediction 
that corrections are 0{e~^^). 

In conclusion our data for the g = 20 Potts model are in excellent agree- 
ment with Finite Size Scaling predictions. In particular the data for CV{l3t) 
agree with FSS for all lattices considered. 

This indicates that the disagreement between FSS and data for the 2- 
d g = 10 Potts model were small system effects. It means however than 
the asymptotic behavior predicted by Binder, and later proven by Borgs 
and Kotecky, only sets in for very large lattices. The lattice size must fulfil 
the conditions L >> ^, L"'"^ >> l/A°'^ (If Widom's relation [0 holds 
this condition is equivalent to the first one), and L'^ » CojiEo — EaY, 
L'^ » Cd/{Eo — EdY, where >> means five to ten times larger. It is 
unfortunate that for such large systems, the transition is already blatantly 
first order. 
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L 


M sweeps 


bins 


1 fi 


^ 7 


s 

O 


16 


10.2 


20 


16 


90.5 


180 






90 


18 


10.0 


20 


18 


10.0 


20 


18 


30.0 


60 


20 


60.0 


120 


24 


7.7 


15 


24 


10.2 


20 


24 


50.5 


100 


24 


20.2 


40 


24 


40.5 


80 


30 


10.2 


20 


30 


10.0 


20 


30 


50.5 


100 


32 


40.5 


80 


34 


7.2 


14 


34 


29.7 


59 


38 


40. 


80 


38 


10. 


20 


38 


10. 


20 



Table 2: ^ = 20 2-d Potts model: Simulated lattices. Lattice sizes, number 
of mega sweeps performed and number of bins used to perform the analysis. 
The first bin may contain a different number of sweeps that the others. It 
has always been discarded for thermalization. 



9 



L 








CV{f3t)/L'' 


16 


1.03647 


(15) 


-0.59302 


(15) 


1.07660 


(6) 


0.2066 


(12) 


18 


1.03598 


(17) 


-0.58489 


(17) 


1.06186 


(5) 


0.2010 


(16) 


20 


1.03582 


(19) 


-0.57878 


(18) 


1.05080 


(5) 


0.1988 


(22) 


24 


1.03482 


(10) 


-0.57013 


(10) 


1.03602 


(2) 


0.1924 


(17) 


30 


1.03360 


(13) 


-0.56230 


(14) 


1.02352 


(2) 


0.1952 


(34) 


32 


1.03317 


(22) 


-0.56058 


(24) 


1.02080 


(3) 


0.1923 


(58) 


34 


1.03286 


(18) 


-0.55911 


(17) 


1.01851 


(4) 


0.1837 


(55) 


38 


1.03218 


(12) 


-0.55653 


(12) 


1.01492 


(2) 


0.1941 


(56) 



Table 3: g = 20 2-d Potts model: Results for CVmax/L"^, BLmin, UAmin and 
CV{Pt)/L\ 
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Figure 1: CVyyiax/L"^ as a function of for the 2-d q=20 Potts model. 
The three curves are drawn using the central value and the one standard 
estimates for Co- 
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Figure 2: BL^nin as a function of for the 2-d q=20 Potts model. The 
three curves are drawn using the central value and the one standard estimates 
for Co. 
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Figure 3: UAjnin ^ a function of for the 2-d q=20 Potts model. The 
three curves are drawn using the central value and the one standard estimates 
for Co. 
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Figure 4: (3{CVmax) and (3{BLmin) as a function of l/Li^ for the 2-d q=20 
Potts model. The curves are drawn using the central value for Co- The effect 
of the uncertainty on Co would not be visible. 
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Figure 5: CV{l3t)lL'^ as a function of l/L^ for the 2-d q=20 Potts model. 
The three curves are drawn using the central value and the one standard 
estimates for Cg- 
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